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Introduction 

The defect structure of B2 FeAl has been the subject of several experimental [1-4] and more numerous but 
conflicting theoretical studies [5-12]. Some support the contention that FeAl is a triple defect forming compound. For 
Fe-rich alloys, the Al-sublattice is always fully occupied, therefore, when an excess of Fe occurs it will substitute on 
A1 sites with no signiflcant composition dependent vacancy concentration occurring on dther sublattice. But with 
excess Al, vacant sites or constitutional vacancies are formed in the Fe-snblattice. For example, the lattice parameter 
measurements for FeAl as a function of stoichiometry [1-4] show a peak occurring at the stoichiometric Fe-50 at. % Al 
composition, udth a linearly decreasing lattice parameter on mther side of stoichiometry. This behavior is consistent 
with a triple defect forming compound [5]. Further arguments in favor of a triple defect structure are based on the 
large heat of formation of FeAl [5,6] and early theoretical work by Nenmaim et aL [7,8] using a pair-wise interacUons 
model for the study of point defects also favors the formation of triple defects in FeAL 

Other studies support the notion that FeAl is an antistmctnre or snbstitntional defect structure compound, where 
the element in excess locates on the lattice site of the minority dement. Work by Kim [9], Weber et aL[10] and Ho 
and Dodd [4] indicate that structural vacancies in the FeAl system are extremdy unlikdy and that the presence of 
thermal vacancies alone can satisfactorily account for the composition and temperature dependence of the vacancy 
concentration in FeAl alloys. Furthermore, the results of Hosoda [11] and the flrst-prindples calculations by Fn et 
al. [12,13] indicate that a triple defect structure modd does not adeqnatdy describe the structure of FeAl and that 
antisite ddects are preferred at the transitional-metal rites for Al-rich aRoys. These results, suggesting that FeAl is a 
substitutional defect forming compound, are more conristent with the observed slip vector in FeAl alloys. 

The absence of substitutional antisite defects at the transition metal sites in confirmed triple defect compounds 
such as NLAl [5,6,12,14,15] indicates that Al atoms prefer to avoid other Al atoms at the nearest neighbor distance. 
This is consistent with the usual (100) slip vector in NlAl [14] since a |(111) partial slip vector would result in nearest 
ndghbor contact between the same type of atoms. FeAl on the other hand deforms by a (111) slip vector [16], at least 
at room temperature, which is much more consistent udth an antirite defect structure. 

Consequently, a strong case can be made for dther a triple defect structure or an antisite defect structure in FeAL 
Yet, alm ost all of these studies fril to point out that no matter what the defect structure is on the Al-rich ride of 
stoichiometry, it is not stable over a significant compositional range as any cursory look at the FeAl phase diagram 
[1,6] will indicate. The B2 structure of the FeAl compound exists to less than a two percent Al-rich deviation in 
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l^ble 1: LMTO results, ECT parameters for bcc-based Fe and A1 and BFS parameters for FeAl 


stoichiometry. This would indicate that neither defect structure could be that stable in the Al*rich end of the FeAl 
system. 

Therefore, the objective of this paper is to provide a consistent description of the defect structure of FeAl by 
modelling the energetics of this system at zero temperature using the BFS method [15,17]. Several basic issues are 
covered in this study including determination of the B2 phase field, determination of the defect structure on both sides 
of stoichiometry, and modelling of the ^ttice parameter as a function of composition. Furthermore, the energetics of 
the BFS model provide a rational for the limited structural stability of the B2 compound for Al-rich compositions. 

The BFS Method 

In this work we follow the same procedure used in a previous appUcation of BFS to the analysis of the defect 
structure of NiAl alloys [15]. For the sake of brewty and in order to allow for a more detailed discussion of the 
results, we refer the reader to pretdons papers that provide a thorough description of the method and its appCcation 
to specific problems [15,17]. In what follows we provide a brief description of the BFS method sufficient to understand 
the discussion of the results that follows. 

The BFS method is based on the idea that the energy of formation of an alloy is the superposition of individual 
contributions of non-equivalent atoms in the alloy [17]: £,• = cf + 9i{£f — ef*). The first term is the strain energy 
£{, computed with equivalent crystal theory (ECT) [18], that accounts for the actual geometrical distribution of the 
atoms surrounding atom t, computed as if all its nrighbors were of the same atomic spedes, and a chemical energy 
£^ — £^ , which takes into account the fact that some of the neighbors of atom i may be of a different chemical spedes. 
For ef we interpret the chemical composition as a defect of an otherwise pure crystaL We represent this defect by 
‘perturbing* the electronic density in the overlap region between dissimilar atoms and locating them at equffibrium 
lattice rites of atom L To ffee the chemical energy of structural defect energy, which should only be induded in the 
strmn energy, we reference cf to a similar contribution where no such perturbation is induded (£p°). The coupling 
function g{, which ensures the correct asymptotic behavior of the chemical energy, is defined as gi = e~*f, where af is 
a solution of sf = —Eq [l — (1 + of )eip(-of )] [19], and where is the cohesive energy for atom £. In the context of 
BFS, the terms ‘strain’ and ‘chemical’ represent quite different effects than the usually assigned meanings. The BFS 
strain energy is related to the usual strain only in that the atomic locations are those found in the actual aBoy: the 
BFS strain energy of a given atom is then the actual strain that it would have in a monatomic crystal of the same 
spedes of the reference atom. Likewise, the BFS chemical contribution is related to the usual chemical energy in that 
the actual chemical composition of the alloy is taken into account, but with the ndghboring atoms located in ideal 
atomic rites: the BFS chemical energy of a ^ven atom is then the actual chemical energy in an ordered environment 
with a lattice spacing characteristic of the equiUbrium lattice of the reference atom. The parameters needed by the 
BFS method, (induding the ECT parameters) were calculated using the Linear Muffin-Tin Orlntals (LMTO) method 
[20] in the Atomic Sphere Approximation (ASA) and are listed in l^ble 1. For brevity, we refer the reader to Ref. 18 
for a detailed discussion on the parameterization of the BFS method. 
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Results and Discussion 

The BFS calculation was performed on a 72 atom cell, allowing for isotropic atomic relaxation induced by the 
presence of vacancies and substitutional atoms. Instead of seafching_tor an absolute energy minlmum.for a ^ven 
composition, we chose to construct a large number of ‘candidate', high symmetry, distributions to obtain information 
on the energetics of the system close to the ground state. This is done with the intention of identifying possible 
metastable states, in addition to the ones corresponding to ideal thermal equilibrium. For example, if configurations 
Cl and Ca both have negative heats of formation, with AiTc, < AJc^, then configuration Cj might still appear in 
the actual alloy with a certain probability, which is lower than that for Ci. If the set of configurations sampled is 
sufficiently large and the structures are chosen respecting the symmetries that characterize the system, one would 
expect to find the ground state, or states sufficiently close to it, for each composition. In order to concentrate on the 
focus of this paper, we leave the details on this computer simulation approach for a forthcoming publication. Examples 
of the atomic distributions for many of the candidate states are illustrated in Fig. 1, and include combinations of Fe 
and A1 antisite defects as well as Fe and A1 vacancies. 

Fig. 1 also displays the heat of formation as a function of A1 concentration for the set of ‘candidate’ configurations 
examined. The results of the simulation are summarized as follows: 

1) The BFS results predict the range of composition for the B2 phase field in excellent agreement with experiment. 
No negative energy states were found beyond 52 at. % Al. There seems to be no limit on the Fe*rich side where the 
heat of formation first decreases in magnitude and then increases again for higher Fe contents, possibly hinting at the 
eristence of the FesAl ordered (DO 3 ) structure. For those compositions consisting of stoichiometric FeAl or Fe-rich 
alloys, the lowest energy configurations were those composed of anitsite defect structures. 

2) For Al'rich alloys, the antisite defect (Al in Fe sites) was found to have the lowest energy (configuration 
(b)). A few other configurations with only slighlty smaller values for the heat of formation were found ((c)-(d)). 
These configurations consisted of arrangements of structural vacancies with only a slightly smaUer energy than for 
the antisite defect structure. All Al-rich structures exhibited much smaller energies than configurations on the Fe-rich 
side of stoichiometry. The sudden termination of the B2 field at xai ~ 52 at. % can be explained in terms of the 
BFS contributions to the energy of formation as described in Eq. (1), and shown in Fig. 2. In this figure, the total 
BFS strain and chemical ener^es per atom are displayed for the lowest energy states for each compoation. The total 
energy of formation is also shown. The BFS strain contribution shows a very slight monotonlc increase with Al content 
throu^ont the whole range of the B2 field, due probably to the greater abundance of large Al atoms. In contrast, 
the BFS chemical energy shows a discontinuous jump for xju > 50 at. %. In this region, substitutional Al atoms 
occupying Fe sites eliminate chemically favorable (a negative contribution from tP — eP> in Eq. ( 1 )) Fe-Al bonds, while 
at the same time increaring the (always positive) BFS strain due to their larger size. The increasing strain reduces the 
magnitude of the coupling function g (which links the strain and chemical components) thus emphasizing the loss in 
negative chemical energy. These combined effects drive the total energy of formation up, quickly saturating the ability 
to absorb additional Al substitutional atoms. On the Fe-rich side, the opposite is true: substitutional Fe antisite atoms 
eliminate favorable Fe-Al bonds as well, but the loss in negative chemical energy is offset by the decreasing BFS strain 
and the subsequent increase in the magnitude of the coupling function. 

3) Being a T=0 K calculation the ground state for stoichiometric FeAl naturally corresponds to a perfectly ordered 
B2 structure. The only other configurations with (much smaller) negative heat of formation corresponds to an antisite 
defect structure. The triple defect structure (two vacant Al sites plus an Al antisite atom) is even more unlikriy, with 
a neariy positive heat of formation. The small value of the heat of formation for this type of defect, relative to the 
corresponding one for the ordered B2 structure or even the antisite structure, hints to the unlikelihood of finding this 
defect in the actual alloy. Consequently, these results are conristent with other theoretical predictions such as Fh et al. 
[12,13] which conclude that the triple defect structure is energetically unfavorable, but our results also indicate that 
the triple defect structure is not energetically impossible. 

4) For Fe-rich FeAl the substitutional Fe-antisite atom is the main defect type, in agreement with all prerious 
studies. Several configurations of Fe-antisite defects, very close in energy, exist for each composition. The difference 
between these configurations resides in the relative position of the Fe antisite atoms. As shown in Fig. 3, where diagrams 
showing the configurations studied for xju = 44.44%, slight gains in energy are obtained as the additional Fe atoms 
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accomodate themselves so as to maximize the correlation between them. In other words, the energy slightly lowers as 
the Fe atoms agglomerate. With even farther deviations from stoichiometry the gap between different configurations 
widens with decreasing A1 content indicating the preference for a particular ordering pattern for the Fe-antisite atoms 
(as the one found for Fes Al, which correspofids~Krthe"D03 ordered structure). 

5) A decrease in lattice parameter with Fe content is also found in this calculation, in agreement with the the 
trend found experimentally. While there is very good agreement for Al-rich alloys, the rate of decrease of the lattice 
parameter noth increasing Fe content is higher than the one observed experimentally. Fig. 4 compares the experimental 
values of Refs. 1-4 with the ground-state values for the theoretical lattice parameter obtained with BFS. The difference 
between the measured and computed values of the lattice parameter as a function of composition can be explained in 
terms of several arguments. Including the fact that this calculation was done at T=0 K and that no local relaxation 
effects (i.e. small displacements of individual atoms around the defect) were included. In addition, consistent with the 
previous discusrion, the clustering of several structures (as shown in Fig. 3) with slightly different order patterns close 
to the ground state suggest that a fair comparison to experimental results should include the effect of these structures 
which individually induce an increase in the lattice parameter. Finally, the actual lattice parameter measurements may 
be influenced by the thermal vacancies that are quenched into FeAl alloys in large numbers [4,11-13] but are not taken 
into account in our calculations. These factors, when accounted for, could only improve the quality of the agreement 
between experiment and the present theoretical analysis. 

Conclusions 

A good understanding of the factors influencing the atomic structure of FeAl or any other ordered system is nec- 
essary if further conclusions regarding physical and mechanical properties are to be derived from any given calculation. 
Our calculation for FeAl seems compatible with most of the evidence gathered through the years, agreemg on the 
essential facts (correct prediction of the B2 phase field, antisite substitutional defects for all comporitions, decrease of 
the lattice parameter with increasing Fe content) as well as ^ving some indication on the probability of other defect 
structures. 
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